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Asymmetry of the natural line profile for the hydrogen atom. 
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The asymmetry of the natural line profile for transitions in hydrogen-like atoms 
< >— ■ ^ ' is evaluated within a QED framework. For the Lyman-alpha Is — 2p absorption 

transition in neutral hydrogen this asymmetry results in an additional energy shift 
of 2.929856 Hz. For the 2si/2 — 2p 3 / 2 transition it amounts to -1.512674 Hz. As a 
new feature this correction turns out to be process dependent. The quoted numbers 
refer to the Compton-scattering process. 
PACS numbers: 31.30.Jv, 12.20.Ds, 0620Jr., 31.15.-p 



The problem of the natural line profile in atomic physics was considered first in terms of quantum 
mechanics by Weisskopf and Wigner JJJ . Within the framework of modern QED it was first formu- 
lated for one-electron atoms by Low (g| . In [g| the occurance of the Lorentz profile in the resonance 
approximation was described and the non-resonant corrections were estimated. Later the line profile 
QED theory has been modified also for many-electron atoms || (see also ]4 1 ) and has been applied to 
the theory of overlapping resonances in two-electron highly charged ions [ I . Within this formalism 
corrections to the energy levels, e.g. reference state QED corrections, have been evaluated ||. 
fi , One of the important consequences of the line profile theory is the occurance of non-resonant 

corrections S. These corrections indicate the limit up to which the concept of the energy of an 
^-j, excited atomic state has a physical meaning - that is the resonance approximation. The exact 
t£^ • theoretical value for the energy of an excited state defined, e.g., by the Green function pole, can be 
^~-f, compared directly with measurable quantities only within the resonance approximation, for which the 
line profile is described by the two parameters: energy E and width T. Beyond this approximation 
t-H \ the evaluation of E and F should be replaced by the evaluation of the line profile for the particular 
^ • process. If the distortion of the Lorentz profile is still small one can formally consider the non- 
resonant correction as a correction to the energy shift. Unlike all other energy corrections, this 
.— ». ■ correction depends on the particular process under consideration which has been employed for the 
r~^. , measurement of the energy difference. The non-resonant (NR) corrections were considered for H-like 
(~^) ■ ions of phosphorus (Z = 15) and uranium (Z = 92) in M, ||. While for uranium the NR correction 
was found to be negligible, its value was comparable with the experimental inaccuracy in the case 
of phosphorus. 
c/3 , In this paper we demonstrate that the concept of a transition energy depends on the measurement 

process. For this reason we evaluate the NR corrections for the neutral hydrogen atom. We consider 
the process of the resonance Compton scattering as a standard procedure for the determination of 
J^h * 1 the energy levels. For this process the parametric estimate of the NR correction can be expressed 
O '_ as H (in rclativistic units). 



Cma 2 (aZf (1) 



5_J ■ where C is some numerical factor, a is the fine structure constant, Z is the nuclear charge number. 

The recent QED calculations for low-Z H-like atoms incorporate corrections of the order of 
ma 2 (aZ) 5 B - [Q; corrections of order ma 2 (aZ) 6 \n (aZ) are also included pi, [H. For \ow-Z 
accurate Lamb-shift calculations to all orders in Za have been performed recently ]15| ], ua . 

Thus in principle the next order corrections to the energy levels should include NR corrections 
and will depend on the process of measurement under consideration. However the numerical factor 
C in Eq. (1) appears to be quite small: 10~ 3 for the Lyman-alpha transition (see below). There 
arc no direct measurements of the Lyman-a transition frequency with an accuracy required in 
order to observe the correction Eq. (|l|). Modern experimental techniques employed in Lamb-shift 
measurements are based on two-photon resonances, e.g. for the transition 2s -Is [17] - [19]. Although 



the theoretical evaluation of the NR corrections to the two-photon resonances is more involved, we 
can state that the corresponding NR corrections will be of the same order of magnitude as the one 
considered here. 

Consider the process of photon scattering on a one-electron atom. This process is described by 
the two Feynman diagrams of Fig. 1. Resonance scattering implies that the frequency of the initial 
photon us is close to the energy difference u> — Ea' — E B where A' is some excited atomic state, B 
is the initial state. Within the resonance approximation we retain only the term n = A 1 in the sum 
over intermediate states in the amplitude, corresponding to Fig. la. The amplitude, corresponding 
to Fig. lb, has a non-resonant character. Separating out the resonant term we express the amplitude 
corresponding to Fig. 1 in the form 

U AB (ujm\; ujm\) = e J — 

^-v (lA UJ ,j, m ,y)An{'jA u ,j m \) n B 

iri, E n - E B - lu 
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where 7 are Dirac matrices, AujmX is the vector potential of the electromagnetic field (photon wave 
function) and E n are the one-electron energies. 

The second and third term in Eq. (2) represent the non-resonant corrections to the scattering 
amplitude. 

The Lorentz line profile arises when we sum up all the electron self-energy insertions in the internal 
electron line in Fig. la within the resonance approximation J2J. After the summation of the arising 
geometric progression one finds 
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where S^(o; + Ea) is the renormalized electron self-energy operator. 

In the resonance approximation we can substitute lu + Ea = Ea 1 in the denominator of Eq. (3). 
The real part of the matrix element Yjr{lu + Ea))a'A' yields the lowest-order contribution to the 
Lamb shift while the imaginary part which is finite and not subject to renormalization indicates the 
total radiative (single-quantum) width of the level A': 

t R {u + E A )) A , A ,=L^-^T A ,. (4) 

An additional contribution to the lowest-order Lamb shift L\T originates from the vacuum polar- 
ization graph. However, this graph gives no contribution to the width [H. 

Taking the square modulus of the amplitude Eq. (2), integrating over the directions of the 
absorbed and emitted photons and summing over polarizations we obtain the Lorentz profile for the 
absorption probability 

jw/ ^ 1 ^A'Adco 

dw ^ = ^^ A , A + L A --.r + i/4r A ,- (5) 
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Here dW(uj) is the probability for photon absorption in the frequency interval [u),u + cLo], 
Ea> — Ea and Ta'a is the partial width of the level A' , connected with the transition A' — * A. The 
inclusion of the Lamb shift La> corresponding to the initial state A into the Lorentz profile (5) can 
be accomplished, if necessary, by the methods developed in []3j, |4|]. The line profile for the emission 
process A' — > A is described again by Eq. (5). 

In the non-resonant terms in Eq. (2) we can substitute 10 — Ea 1 — Eb- Then we arrive at the 
expression (we omitted the Lamb shift La>, which is not essential for our purposes) 
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The differential cross section of the process is 

da AB {ujmX; u'j'm'X') = 2n\U < £ B (wjm\; Lu'j'm'X')\ 2 5{E A -E b +lu' - u)duduj' (7) 

The onc-clcctron states A,A',B depend on the usual set of quantum numbers A = n A l A j A m A , 
where ua is the principal quantum number, I a is the orbital angular momentum quantum number 
which determines the parity, j A , ttla are the total angular momentum and its projection. Integration 
over u) 1 , summation over j'm'X' , ttia and averaging over to, m B yields : 
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Substituting Eq. (6) in Eq. (8) and omitting the terms, quadratically dependent on the non- 
resonant contributions, we obtain 
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The first term in Eq. (9) leads to the usual Lorentz line profile for the process under consideration 

1 T AA ,W BA ,{jX) 
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where W B a> is the transition probability B — > A' 

WbA '^2~ b -TT E \h A *«imx)BA>\ 2 

mmBtn A f 

and Yaa' is the partial width of the level A' , connected with the transition A' — > A 

T AAI =^W AA ,{j'X'). (12) 

It is assumed that the level B is stable (T B = 0) or mctastablc (T B <C T^/). 

The term cr AB that represents the interference between the resonant and non-resonant contribu- 
tions to the amplitude is determined by 
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where Wbb ; „a'C?A) and W nB -AA'{j^) are the "mixed" transition probabilities 
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W nB ;AA> C?A) = — — J^ (7^jmA)An(7-^*wjmA)BA' , (15) 
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and r^sjcn is the "mixed" partial width 

I^B-CD - ^ WUfljCB (/A') . (16) 

j'X 

It is important to emphasize that in the sum over n in Eq. (13) only the states n with the same 
symmetry (i.e. with the same Ja'^A') survive. Thus we keep the same averaging factor (2j A > + 1) _1 
in Eqs. (14), (15). Note also that for the fixed transition B — > A 1 the type of the absorbed photon 
(JX) will be also fixed. 

We assume that the standard way of measuring the resonance frequency is employed which is 
connected with the determination of the maximum in the probability distribution for the given 
process. In the pure resonance case the maximum condition 

^H =Q (17) 

corresponds to the resonance frequency value w max = Ea> —Eb- If we take into account the correction 
a AB (oj), the result will be different 



with 



and 



A (a W M+(T W H ) = o (18) 



Wmax = E A > - E B + S (19) 
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Thus the value of u) max for the probability distribution in the photon scattering process cannot 
be compared directly with the energy difference between the two levels. The process-dependent 
non-resonant correction S should be taken into account. This result holds also if we include any 
QED corrections in E A i , Eb ■ The order of magnitude of the correction S follows from the standard 
estimates for the allowed transition probabilities (also the "mixed" ones) ma(aZ) A and the transition 
energies m(aZ) 2 : 

[ma(aZ) 4 } 2 2 6 

= ma (aZ) (21) 



m{aZ)' 1 

These orders of magnitude are the same whether or not the level E A i is removed from Eb by a 
fine (hyperfine) structure splitting B. 

Consider now the neutral hydrogen atom in the non-relativistic approximation. We put B = 
A = Is, A' = 2p. Then the line profile will correspond to the Lyman-a transition. In this case 
Ta> = r ' AA i = Wba 1 and the matrix elements of the electron-photon interaction are 

{l-Kjj m \)BA' = {lA u j m \)A'B ~* (^lm)lsls , (22) 

where 



U lm = ^\ l /^ 3 /2 rY lm (23) 



and u>o = E 2p — E\ s . In Eq. (23) and below we use atomic units. 
Then the first term of the correction 8 can be written as 



dr ir2 VL (ri) U lm (n )G Ba , (n ; r 2 ) t/^, (r 2 ) Vi s (r 2 ) (24) 
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where Gs., (ri;r 2 ) is the non-relativistic "modified" Coulomb Green function and ipi s is the 
Schrodinger wave function. For the angular integration we use the partial wave expansion for Ge a , 

G EA/ (ri;r 2 )=J2G l E A ,(n;r 2 )Y l * M (a 1 )Y lM (Q 2 ), (25) 

IM 

where Yjm(^) are the spherical harmonics. The angular integration and the summation over the 
angular momentum projections yields : 

-i 6 /*oo />oo 

6l = 2WJ J ^i^r?r 2 3 VL(r-i)G^,(r- i; r 2 )^ 8 (r 2 ). (26) 

We evaluate Eq. (26) with the Sturmian expansion for G E , (j"i; r 2 ) that we take from [p0| : 
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where i?„; are the radial Coulomb wave functions. Insertion of Eq. (27) in Eq. (26) and integration 
over r\ , r 2 for Ea> = 2p results in 

* - T (0 16 (± (m ^2), ^< 2 - *»• * 4: 2 ' 3 > + 7 ' 2 ) ■ < 28 > 

where 2 iq denotes a hypergeometric function. The expansion in Eq. (28) converges very fast and 
for m = 10 it gives an error less than 10~ 6 . We obtain 

81 = 2.127209 • 10~ 3 a 6 = 2.1168998 Hz . (29) 

The second term of the correction 8 can be written again as Eq. (26) but with the "normal" 
Coulomb Green function Ge a +e b -e a , (^i;r 2 ) = G_ 7 / 8 (ri; r 2 ). In this case it is convenient to use 
for the radial Green function G^(ri;r 2 ) the representation |2(J 



G l E (r i; r 2 ) = - V -^^R ml (2 ri /v)R ml {2r 2 /v) (30) 

v * — ' m — v 

m=2 

where v — Z/y / —2E = 2j\fl. Now we obtain 



7 00 

(in + 

m — v 



h = ^j^t ( ~ ) V (m + lMm " 1} 2 ^ 2 (2 - m, 5; 4; 2/(* + 1)) . (31) 
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Retaining only six terms of the expansion (31) yields an accuracy of 10 6 . Then 

S 2 = 0.821625 ■ 10~ 3 a 6 = 0.81764337 Hz (32) 

and finally 

6 u P L =h+$2 = 2.929856 Hz . (33) 

It should be mentioned, that the Lyman-a resonance consists of two peaks in the scattering exper- 
iment, corresponding to the two fine-structure components. In the non-relativistic approximation 
the distortion of these two peaks is equal and defined by Eq. (33). 

We made also an analogous calculations for the transition 2s\/2 —> 2p 3 / 2 ~^ lsi/2- In this case 
B = 2si/2, A = lsi/2, A' = 2p 3/2 . The result is 



S 2 S P ii 2) = - I I ) « 6 = -1-522439 ■ 10~ 3 a 6 = -1. 512674Hz. (34) 
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Concluding, we can state that the NR corrections are comparable with some QED corrections 
to the Lamb shift recently included in the consideration [[L2|. The level of accuracy of modern 
experiments with neutral hydrogen for the two-photon transition 2s — Is is estimated in total to be 
46 Hz [jl9|| and is approaching the magnitude of these non-resonant corrections. The derivation of the 
two-photon NR corrections requires consideration of two-photon resonance scattering on a hydrogen 
atom. The corresponding expression for the amplitude will contain two energy denominators unlike 
as in Eq. (2). However, only one of them will become resonant while the other one appears as a 
nonresonant factor and will not change the order of magnitude of the total amplitude. Accordingly, 
the interference term (13) defining the resonance shift will lead to values of the same order of 
magnitude as the corrections (33) and (34), respectively. 
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FIG. 1. Feynman graph, corresponding to the photon scattering on an atomic electron. The solid double 
line denotes the bound electron in the atom. The wavy lines denote the absorbed and emitted photons. 
The initial, intermediate and final states of an electron are denoted as B, n, A. The initial and final photon 
states are uijrnX and ui'j'm'X' where w is the frequency, jm denote the photon angular momentum and its 
projection, the quantum number A determines the parity of the photon state. 



